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INTRODUCTION 


Accurate theoretical descriptions (refs. 1 to 3) of high-energy heavy-ion scat- 
tering phenomena require knowledge of nuclear single-particle density distributions 
(i.e., matter densities). Experimentally determined nuclear density distributions ^ 
however, are limited to charge and magnetization densities, obtained from electron 
scattering experiments. Extensive tables of nuclear charge distribution parameters 
exist (refs. 4 and 5). In practice, the charge densities are usually directly sub- 
stituted for the matter densities where required in the analyses. The resultant 
theoretical cross sections, however, are typically overestimates (refs. 3, 6, and 7). 
As discussed in reference 3, these overestimates are a result of the increased sur- 
face diffuseness of the nuclear charge distributions over that of the matter distri- 
butions because of the finite spatial extent of the proton charge. In references 2 
and 3, a method for extracting matter density distributions from Woods -Saxon charge 
densities was presented in conjunction with a generalized optical -model reaction 
theory for heavy-ion scattering. Although Woods -Saxon distributions are adequate for 
many nuclei, a better representation for the charge distribution of light and medium- 
weight nuclei (three to eight protons) is the harmonic well distribution (refs. 4, 5, 
and 8) based upon the nuclear shell model (ref. 9) . In this paper, the general 
method of references 2 and 3 is applied to obtain matter densities from harmonic well 
charge densities. Unlike the Woods-Saxon distribution (used in refs. 2 and 3), which 
requires some analytic approximations to obtain the resultant matter density expres- 
sion, the harmonic well distribution is exactly solvable. 

As developed in references 2, 3, and 10, the heavy-ion reaction theory ignores 
the Pauli exclusion principle (Pauli correlation effects) by using simple products of 
wave functions rather than properly antisymmetric ones for calculating nuclear matrix 
elements. Although unimportant in estimating total and absorption cross sections 
for high-energy heavy-ion collisions (ref. 2) , these Pauli correlation effects become 
significant when determining projectile-nucleus abrasion cross sections for small 
residual fragment masses (i.e., when there is large overlap between projectile and 
target nuclear volximes). In this paper, a method for incorporating these effects 
into the heavy-ion reaction model is presented. Finally, cross sections for selected 
nuclei are predicted with the improved optical potential model. 


HARMONIC WELL MATTER DENSITIES 

For relatively light nuclei with incomplete Ip shells, the harmonic well charge 
density (refs- 4, 5, and 8) is given by 


p(r) 



( 1 ) 


where is a normalization constant, r is the radial coordinate, the parameter a 

is some simple function of Z, the nucleus proton number. 


a = f (Z) 


( 2 ) 



and a is the •'oscillator parameter" given by 


a 


(mo) ) 
osc 


" 1/2 


( 3 ) 


where m is the nucleon mass, 938 MeV/c^. In equation (3), ^osc 
equidistant-energy interval between successive oscillator levels (ref. 8). Rather 
than determining a from some specific functional dependence on Z, it is some- 
times treated as a free parameter. Equation (1) is then referred to as a "modified" 
harmonic well charge density. Specific values for a and a, for the nuclei of 
interest, are tabulated in references 4 and 5. Values of the normalization constant 
pQ are obtained by requiring that 

X oo 

r^ p(r) dr = 1 (4) 

The matter density is related to the harmonic well charge density through the 
expression 


p (r) « / p (r*) p (r+r*) d^r’ 
c J ^m 


(5) 


where p«(r) is the nucleus charge density, p«(^) is the proton charge density, 

and is desired nuclear matter density. Equation (5) implies the assump- 

tion that the neutron and proton number densities within the nucleus are identical. 
This is reasonable since they differ only because of Coulomb repulsion between the 
protons which is a small effect in lighter nuclei. Taking the Fourier transform of 
equation (5) and using the convolution theorem yield a simple product of form 
factors: 


F (q) = F (q) F (q) 
c pm 


where 


F(q) 


= r 

^ Jn 


r(sin qr) p(r) dr 


and q is the magnitude of the momentum transfer. 

From equations (1) and (7), the nuclear charge form factor is 


„ , , 3/2 3 (^ . 3a \ -q 

Fc(q) = ^ ~ 4 " ) 


2 2 
a 


( 6 ) 


(7) 


( 8 ) 
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The usual form (refs. 4 and 8) for the proton charge density is a Gaussian function: 


Pp(?) 



(9) 


where r » 0.87 fm (ref. 11) is the proton root -me an -square charge radius. From 
equation^ (7) / the proton charge form factor is 


Fp(q) = — 


2 2 

Hi 


Hence, equations (6), (8), and (10) yield a matter density form factor of 


( 10 ) 


/ 2 2 

„ , , 3/2 3 / , . 3a aq a 

F(q)=p7l' ^ 

m *^o \ 2 4 


2 2 

exp(-q s ) 


where 


2 

s 


2 

a 


4 


2 

r 

P 
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( 11 ) 


( 12 ) 


The matter density is obtained by utilizing the inverse Fourier transform of the 
matter density form factor: 





q(sin qr) F (q) dq 
m 


(13) 


Substituting equation (11) 


into equation (13) yields 





^ 2 2 2 
3aa oa r 

2 4 

8s 16s 


2 

-r 

exp — - 
4 s 


(14) 


Table I lists the values for a and a, taken from reference 5, which are substi- 
tuted into equations (1) and (14) to determine the charge and matter density distri- 
butions. These results are plotted in figures 1 through 6 for six nuclei. 
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PAULI CORRELATION EFFECTS 


Since nucleons have half-integer intrinsic spins (i-e-, they are fermions), they 
must obey the Pauli exclusion principle. This principle requires the overall system 
wave function to be antisymmetric under the exchange of any two nucleons in the 
system. For a system composed of only two nucleons, the properly antisymmetric wave 
function is 


(HI, 2) = (1) <|)-(2) - ^ 12 ) (|)-(1)] (15) 

Up Up 

where and (t)^ are the wave functions for the nucleons. The product, 

4)^(1) 4)p(2) 

is usually called the direct term. Subtracted from the direct term is a product 
representing the exchange of coordinates of the two nucleons. This product. 


U p 

is sometimes called the exchange term. 

Previous versions of the optical-model reaction and fragmentation theories 
(refs. 2, 3, and 10) neglected antisymmetry and used only product wave functions (the 
direct term in eq. (15)). For determinations of heavy-ion total and absorption cross 
sections, neglect of antisymmetry is not an extreme assumption since the largest 
contributions to these cross sections come from peripheral collisions where there is 
little or no overlap between colliding nuclei (see "Results”). Pauli correlation 
effects, however, are significant when predicting projectile abrasion cross sections, 
particularly for small residual mass fragments (when there is significant overlap 
between the volumes of the colliding nuclei) . 

In references 2, 3, and 7, the optical potential operator is 


Vt- S 

0C3 


0C3 


(16) 


where t^. is the two-body transition amplitude between the a~constituent of the 
target and the j-constituent of the projectile. The expression in equation (16) is 
general in that it was derived independently of any specific assumptions regarding 
nuclear wave functions. In references 2 and 3, it was shown that for a simple prod- 
uct of wave functions, equation (16) yields for the optical potential 


W(x 


) = ApA^ y' y* Pp(x+y+t,j,) t(e,y) 


(17) 
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where 


Ap,A,p nuclear mass numbers of the projectile and target 

e two-nucleon kinetic energy in their center of mass frame, GeV 

X relative position vector of the projectile, fm 

y two-nucleon relative position vector, fm 

% collection of constituent relative coordinates for target, fm 

and the average two-nucleon transition amplitude is 

= a4“ ^ 

*P*T aj 


Equation (17) is a valid representation for the direct term, but does not include any 
Pauli correlation effects. In the following sections, equation (17) is generalized 
to include the exchange (correlation) effects. 


Second Quantization Notation 

In terms of second quantization notation (ref. 12), a two-body operator can be 
expressed as 


G 


2 


ik Jim 


(19) 


where by is an annihilation operator and a creation operator for a nucleon in 

the single-particle state y = i,k,jl,m) . The matrix element in equation (19) is 
given by 


(ik|g|jlm) = J* J* d 


3-»- 






(b|^{x2) g(x^ 




( 20 ) 


In this notation, a two-particle state is written as 




( 21 ) 


where |o^ is the vacuum state. Thus, in terms of the two-particle state, the two- 
body operator is 


<ik|G|*n> = I E ^ 2 ^ (i'k'|gU'm')<0|b^b^bJ^.b]^.bj^.b^.b^b^|0> (22) 
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The operator commutator relations are 


b.b + b b. = 0 

1 k k 1 


b!b,^ + b,^bt = 0 

1 k k X 


b.b^ + b.^b . =6., 

X k k X xk 


where is the Kronecker delta. Using these relations and the property. 


bJo> = <0|bj = 0 


enables equation (22) to be written 


<xk|G|j^^ = (ik|g|j!jn) - (ik|g|ml) 


If the initial and final states are the same, equation (25) yields 


<xk|Glik> = (ik|g|ik) - (ik|gjki) 


Thus, the matrix elements of two-particle operators, in second quantization, are 
antisymmetric. 

For an A-body state, the state vector is 


[a> = b^...b^...b^|o> 


Hence, for a two-body operator, equations (19) and (27) give 
<aIg|a> = j H (ik|g|jlm)<A|bj^b^b^b^|A> 

-i V 9m 


which reduces to 


<a|g|a> = ^ 2 [(ik]g|ik) - (ik|glki)] 

i,k=1 
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The Optical Potential 


In second quantization notation, the state vector for the ground state 
projectile nucleus is 




where are projectile nucleon state creation operators • 

state of the target nucleus is 


Similarly, the 


IV = Io> 


where g^ are target nucleon state creation operators. Thus 


W(x) =<P„T„|V ^Ip„T„> = <T |<P„ Iv ^1 p„>|t„> 

0 o' opt' 00^ ^ 0'^ o' opt' 0^' 0^ 


opt ' 


yields 


W(x) = Ei: ( Pk| 1 1 aj)<Tp |<Pp 


which reduces to 


Am 

T P 


W(x) 


]C 2 [(aj|t|aj) 
0E=1 j = 1 


- (aj|tlja)] 


The matrix elements are 


and 


(aj|t|aj) =/d\ y'd^x. <!>.($.: 


(ocjltlja) = / d^x / d^x . 4)^(x ) ())t(x.) t . (x ,x.) (t>.(x ) (x_ 

J aj 1 oc OL :} 3 aj a' 3 ^3 a j 


In reference 7, the direct term in equation (34) was evaluated as 

22 (aj|tlaj) = ApA^ r P^(f„) f d^y p (x+y+? ) t(e,y) 

a j J •' 


of the 

(30) 
ground 

(31) 

(32) 

(33) 

(34) 

(35) 

(36) 

(37) 
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which is equation (17). To evaluate the exchange term, in equation (34)/ we rewrite 
equation (36) as 


= f d^x 4>'^(x ) (t>!(x.) t .(X ,x.) ° 

J a J J a a j j aj a j 3 a a j ^ 


( a j 1 1 1 j a) 


a a ' J J 


Rearranging the wave functions in equation (38) gives 

(«j|t|j« ) *.(J ) t (J .J ) I;” ^ 

“'(b(x) (b.(x.) 

J_a a 3 


Combining equations (34), (35), and (39) yields 


f(x) = Idx Idx, (|)(x)(j)(x) <|).(x.) (j).(x.) t .(x /X.)[1 - 

. J CC J 3 a a ^3 3 ^3 3 «3 a 3 


C .(x /X.)] 

«3 a 3 


with the correlation function given by 


) 4>^{x ) 

^ a 3 3 « 

OCl 

V*«> 


In terras of nuclear single-particle densities. 


p(x. ) = (l>|(x. ) <|). (x. ) 

^ X XI XX 


W(x) becomes 


W(x) = 


) = ]C / p^(x ) / d^x. p_(x.) t .(x ,x.) [1 - C ,(x ,x.)] (43) 

.y a^Tay 3 P 3 a 3 U 3 ocj a j 


At this point, it is assumed that depends only upon the relative posi- 

tion of the a- and j-constituents. Then, the Fourier transform of equation (43) is 


W(q) =I^t„.(e/q) n - c (q)] F^(q) Fp(q) 
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where the fact that t^^^j depends only upon the relative position of the a- and 
j -constituents has been used. Ihe form factors F,p($) and Fp(5) are the Fourier 
transforms of the single-particle densities and transform of the 

correlation function. 

Following references 2, 3 , and 7 , the inverse Fourier transform of a single term 
in equation (44) is taken and then summed over all constituents to give 

W(x) = y y Pp(x+y+f_^) t(e,y) [1 - C(y)] (45) 

where an approximation is made by introducing the average correlation function C 
given by 


C 


-JL_ 

A A “ 
P T 


C . 


(46) 


Equation (45) reduces to the previous result, equation (17), if there are no 
correlation effects (i.e., if C = 0). 


Correlation Function Approximation 

Since C^j in equation (41) depends explicitly on unknown nucleon single- 
particle wave functions (t)j^(^j^)/ the correlation function in equation (40) or (45) 
cannot be determined. Therefore, it must be approximated. One very good approxima- 
tion, even for relatively light nuclei, is the infinite nuclear matter approximation 
(ref. 13) where the nucleon single-particle wave functions are assumed to be plane 
waves. The correlation function, from reference 13, is then written as 


C(y) 


2 ^ 

4 


3 


kpy 


(47) 


where j-| is the spherical Bessel function and kp = 1.36 fm"^ is Fermi wave 
number. In this paper, for ease of analysis, the expression in equation (47) is 
replaced by a Gaussian form 


C(y) 



exp(-0. 18y 


(48) 


since 


|3 


V 2 2 

4 4 

= ^ ^ ^ ^ 


( 49 ) 
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and 


2 2 

2 2 y 4 4 

exp(-kp^y^) = 1 - + O' (kp y ) 


( 50 ) 


agree for small values of k^y (i.e./ when correlations are important) 


RESULTS 

Harmonic Well Density Distributions 

Figures 1 through 6 display results for the charge and matter density distribu- 
tions of lighter nuclei (3 < Z < 8) obtained from equations (1) and (14). Values for 
the charge distribution parameters a and «/ taken from reference 5, are listed in 
table I for each of the nuclei. Analytic expressions for the nuclear matter density 
distributions in figures 1 through 6 can be parameterized as 


p^(r) = Aq(B^ + C^r^) expC-D^r^) 


(51) 


where from equation (14), 


A^ = 




= 


8s' 


3g 3aa 
^ ^ 2 " ^ 2 


8s 


= 


oca 


16s 


1 

^0 ” 2 
4s 


(52) 


Values for these matter density parameters are listed by nucleus in table II. 


Heavy-Ion Total and Absorption Cross Sections 
From eikonal scattering theory (ref. 14), the complex phase function is 


X(b) 


1 

2k 



— oo 


U(^,z) 


dz 


(53) 
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where b is projectile impact parameter, k is wave niimber, z is magnitude of the 
projectile position vector in the beam direction. The reduced potential U in terms 
of the optical potential W(^) is 


U(x) = 2mApA^(Ap + A^)""' W(x) 


(54) 


Thus, the phase function, including Pauli correlation effects, is 


x(^) = j ApA^ a(e) [a(e) + i] 1(5) 


(55) 


with 


Kb) 


= [2n B(e)] 


-3/2 


f I 


6T 




/ 


3 ”^ 

d y Pp(b+z+y+^^) 


[1 - C(y)] exp 


-y 

2 B(e) 
(56) 


Values for the nucleon-nucleon scattering parameters a(e), cr(e), and B(e) in 
equations (55) and (56) were taken from the compilations in references 15 and 16 and 
averaged over the projectile and target constituent types as in references 2, 3, 7, 
and 10. When computing cross sections, the correlation^f unction approximation from 
equation (48) was incorporated into equation (56) for C(^). 

In terms of the complex phase function absorption cross section is 

{l - exp[-2 Im ^ 

Table III displays representative absorption cross sections (from eqs. (55) to (57)) 
for carbon projectiles, at several different incident kinetic energies, colliding 
with various target nuclei. Also listed, for comparison, are theoretical predictions 
from reference 2 and available experimental data (refs. 17 through 20). The improved 
agreement between experimental results and the theoretical predictions in this work, 
over the predictions of reference 2, is primarily due to the use of the more precise 
harmonic well densities rather than the approximate Woods-Saxon distributions for 
lighter nuclei. Table IV shows that the effects on the absorption cross sections of 
including Pauli correlation were minimal. In general, the percentage reductions in 
cross sections were largest for lighter nuclear systems (low mass numbers) at 
lower incident kinetic energies. For a given collision pair (Ap, A.p) , the percentage 
reduction decreases as the incident energy increases. For a given energy, the 
percentage reduction also decreases as system mass numbers increase. In all cases, 
however, the reductions in predicted absorption cross sections were less than 
10 percent. 

Table III reveals that the theoretical absorption cross sections from refer- 
ence 2 are consistently smaller than those obtained in this work. This is a manifes- 
tation of the smaller surface thicknesses of the approximate Woods-Saxon densities. 
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used in reference 2, over the actual harmonic well surface thicknesses* The approx- 
imate Woods-Saxon densities used in reference 2 were obtained from more exact 
harmonic well distributions by artificially flattening the central charge densities 
(see ref* 4) at the ej^ense of the longer range tail* The results of this artificial 
flattening for are displayed in figure 7* Note that the 90% to 10% region is 

smaller for the Woods-Saxon distribution than for the harmonic well distribution. 

Figure 8 displays predicted absorption cross sections from this work and refer- 
ence 2 for projectiles, at 2.1 GeV/nucleon, versus target mass number A^* Also 
displayed are ej^erimental results obtained by three e3q>erimental groups (refs. 18, 
19, and 21). Again, the newer theoretical predictions are in better agreement with 
experiment than the predictions given in reference 2* As with the carbon projectile 
results of table III, the cross sections from reference 2 are smaller because of the 
artificial flattening of the oxygen charge density into a Woods-Saxon shape* 

Collision total cross sections are obtained from 


tot 




= 4711 {1 - exp[-Im x(t>)] cos[Re x{h)]} b db 


(58) 


where the phase function, as before, is determined from equations (55) and (56). 

Table V gives results from equation (58) for collisions at two different 

laboratory kinetic energies* Also shown are theoretical predictions from reference 2 
and experimental results from reference 17. The lack of agreement between the pre- 
dictions of this work and the experimental results of reference 17 may be due to the 
large uncertainty in a(e) (around 60 percent at these energies) in Re This 

uncertainty does not affect cr^bs t)ecause a(e) does not appear in equation (57). 
Varying a(e) between its limits of uncertainty (see ref. 2) varies cr^ot 
12c_12c collisions, at 2*1 GeV/nucleon, by 4 percent* This variation is sufficient 
to give agreement with the limited results of reference 17* Improvements to this 
phase of the theory would require either more accurate knowledge of a(e) or more 
experimental data for coirparison, or both. 


Projectile-Ion Abrasion Cross Sections 
From reference 10, the cross section for abrading n projectile nucleons is 


a 

n 



- exp[-A^ a(e) 1(b)] exp[-A^A^ cr(e) 1(b)] b db 


(59) 


where the residual fragment (prefragment) mass number is 


A = 
F 


A 


P 


n 


( 60 ) 


and I(^) is given by equation (56). Results obtained from equation (59) are 

2 0 

displayed in table VI for Ne projectiles at 2* 1 GeV/nucleon colliding with 
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1 2 

targets. Also displayed are abrasion results without Pauli correlations (i.e.^ 

C = 0 in eqs. (56) and (59)). For n < 12, correlation effects are insignificant. 
This results from the peripheral nature of these collisions where there is little 
overlap of the colliding nuclear volumes. As the number of abraded nucleons 
increases, greater volume overlap is required which results in larger correlation 
effects. The Pauli correlation contribution should be greatest when there is 
complete overlap of the colliding nuclear volumes. This situation occurs for 
n = Ap = 20 (i.e., Ap = 0) where, from table VI, the cross section is reduced 
by an order of magnitude from the uncorrelated result. Recall that equa- 

tion (16) is written in terms of free nucleon-nucleon transition amplitudes. In 
reality, the effective force in nuclear matter is much weaker than the free nucleon- 
nucleon force (ref. 22) . This is manifested in the reduced cross sections obtained 
when exclusion principle correlation effects are included. For the results presented 
herein, the carbon target density parameters are those listed in tables I and II. 

For the ^^Ne projectile, the matter density is extracted, using the procedure in 
reference 2, from a Woods -Saxon charge distribution: 


Pc(r) 


Pq 

1 + exp[(r - R)/c] 


(61) 


with nuclear half -density radius R and nuclear surface thickness t given by 
(ref. 23) 


R = 2.740 fm 
t = 2.515 fm 


and diffuseness c obtained from 


c = t/4.4 


(62) 


Until recently, no experimental abrasion data were available for comparison with 
theory. In the past, theoretical abrasion cross sections (which do not include 
important ablation effects) could only be compared with experimental fragmentation 
results (which do include ablation effects) . Hence, experimental verification of a 
particular theoretical abrasion model was not possible. Recently, Stevenson et al. 
(refs. 24 and 25) have provided e3q)erimental abrasion measurements suitable for com- 
parison with theory. Figures 9 and 10 display these experimental results for ^^Ne 
projectiles (at incident kinetic energy of 2.1 GeV/nucleon) colliding with a carbon 
target. 

The e3q>erimental results are presented as relative probabilities for forming a 
particular projectile fragment residual mass Ap by abrading n nucleons from the 
incident projectile nucleus. These displayed results also include only those frag- 
ments with nuclear charge Z less than 10. This minimizes the likelihood of uninter- 
acted beam particles being included in the data. The displayed mass distribution 
does not cut off sharply at Ap = 20 because of the finite resolution («1.5 amu) of 
the detector (ref. 25). To compare theory with these experimental results, the theo- 
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retical abrasion cross sections must be converted into relative probabilities. Let 
the relative probability be 


R.P. 


o 

n 


a* 

abs 


(63) 


where is obtained from equation (59) and 


a' -a - 0.5a 
abs abs 1 


In equation (64), ^abs given by 


abs 


P 


a 

n 


(64) 


(65) 


1 9 

and the term, 0.5a-|, accounts for the Ne fragments which are not included in the 
experimental data. Additionally, for Ap = 19, the relative probability is given by 


R.P. 


0 .5a^ 

a' 

abs 


( 66 ) 


1 9 

since the Ne fragments must be excluded. 

As displayed in figure 9, the agreement between this theory and experiment is 
very good. Note also that the theory agrees better with the experiment when Pauli 
correlation effects are included. This is especially true for the lighter residual 
mass fragments. The slight disagreement between the theoretical predictions (with 
correlations) and experimental results may be indicative of the approximate nature of 
the infinite nuclear matter approximation to the actual correlation function. 
Unfortunately, improvements in this area are hindered by the lack of knowledge of the 
exact nuclear single-particle wave functions. 

Finally, the marked improvement in the abrasion predictions of this work over 
those of the original version of this abrasion theory is displayed in figure 10. 
Recall that the original theory (ref. 10) used approximate Woods-Saxon densities for 
and did not include correlation effects. Surprisingly, both calculations yield 
comparable absorption cross sections (1057 mb for the theory in ref. 10 versus 
1076 mb for this work). 


CONCLUDING REMARKS 

In this work, a previously developed heavy-ion reaction theory, capable of pre- 
dicting total, absorption, and abrasion cross sections for heavy-ion collisions, has 
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been significantly improved in two ways. First/ the accuracy of theoretical predic- 
tions for lighter nuclei { three to eight protons) was improved by developing an exact 
analytic ejqjression for the correct harmonic well matter density distributions rather 
than approximating these by inexact artificially flattened Woods-Saxon distributions. 
Second, the need for Pauli exclusion principle correlation effects to be included 
in the theory has been identified and the theoretical framework extended to include 
them. Since exact analytical methods for calculating correlation functions require 
knowledge of the, as yet, unknown nuclear single-particle wave functions, an infinite 
matter approximation was satisfactorily implemented. These Pauli correlation effects 
were noted to be especially important for predicting those abrasion cross sections 
resulting from collisions involving large spatial overlap of the colliding nuclear 
volumes. Further theoretical improvements are presently being developed. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
March 5, 1982 
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SYMBOLS 



F(q) 

G 

^a'^a 

I(^) 


Im x(^) 


pj'Pj 

q 


Re x(^) 




nuclear mass number 

—3 

defined in equations (52), fm 
oscillator parameter, fm 

average slope parameter of nucleon-nucleon scattering amplitude, fm"^ 

defined in equations (52) 

projectile impact parameter vector, fm 

annihilation operator for Jlth single-particle state 

creation operator for Ith single-particle state 

correlation function between a- and j-constituents 

average correlation function 

-2 

defined in equations (52) , fm 
speed of light, m/sec 

-2 

defined in equations (52) , fm 

two- nucleon kinetic energy in their center of mass frame, GeV 
nuclear form factor 

unspecified two-body operator in equation (19) 

annihilation and creation operators for target a-constituent state 
defined in equation (56) 

imaginary part of eikonal phase shift function 
number of abraded nucleons 

annihilation and creation operators for projectile j-consituent state 
momentum transfer vector, fm"^ 

real part of eikonal phase shift function, dimensionless 
position vector, fm 
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I 


I 


s 

t 

aj 

^opt 

W(^) 

? 

z 

i 


( 



proton root-mean- square charge radius, fm 
defined in equation (12) 

average two-nucleon transition amplitude, MeV 

two-nucleon transition operator for nucleons a and j, MeV 

optical potential operator, MeV 

optical potential (defined in eqs. (17) and (45)), MeV 
relative position vector of projectile (ic = ^ + z) , fm 
two- nucleon relative position vector, fm 
total number of nuclear protons 

position vector of projectile in beam direction, fm 
binomial coefficient 


a harmonic well distribution parameter (see eq. (2)) 

average ratio of real part to imaginary part of nucleon-nucleon scattering 
amplitude 

collection of constituent relative coordinates for target, fm 
p nuclear density, fm"^ 

Pq normalization constant in equation (1), fm“^ 

<7(e) average nucleon-nucleon total cross section, fm^ or nib 

^abs heavy-ion absorption cross section, fm^ or nib 

‘^exp experimental heavy-ion cross section, fm^ or nib 

cro^s section for abrading n nucleons, fm^ or mb 
^tot heavy-ion total cross section, fm^ or mb 

(|) nuclear single-particle wave function (fm)“^/^ 

X(^) eikonal phase shift function 

cp nuclear two-body wave function (fm)"^ 
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Subscripts: 


c 

charge 

F 

prefragment 

m 

matter 

P 

projectile 

P 

proton 

T 

target 


Arrows over symbols indicate vectors. 
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TABLE !•- NUCLEAR CHARGE DISTRIBUTION PARAMETERS 


FROM ELECTRON SCATTERING DATA 
[From ref* 5] 


Nucleus 

— 

a 

a, fm 

7l1 

0.327 

1 .77 

^Be 

.611 

1.791 

^ ^B 

.811 

1.69 

'2c 

1.247 

1.649 


1 .291 

1.729 

'^0 

1.544 

1.833 


TABLE II.- HARMONIC WELL MATTER DENSITY DISTRIBUTION PARAMETERS 


Nucleus 

Aq , fin“2 


Cq, fm"2 

Dq , fm"2 


0.0282 

0.906 

0 . 148 

0.380 

^Be 

.0211 

.829 

.268 

.370 

"b 

.0225 

.739 

.419 

.425 

'2c 

.0190 

.571 

.691 

.452 

'% 

.0156 

.607 

.625 

.402 

'6o 

.0112 

.591 

.636 

.350 
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TABLE III.- ABSORPTION CROSS SECTIONS FOR ^^C PROJECTILES 
COLLIDING WITH VARIOUS TARGET NUCLEI 



^abs' 

mb 





cTexp^ mb 


This work 

Ref. 2 



0.87 GeV/nucleon 


12 

819 

763 

*939 ± 49 


2.1 GeV/nucleon 


1 

237 

246 

*269 ± 14 




^258 ± 21 

12 

839 

781 

^888 ±50 




^826 ± 23 

16 

990 

820 

^1022 ± 25 

64 

1727 

1656 

*^1730 ± 36 

138 

2519 

2447 

*^2600 ± 100 

184 

2924 


■^3000 ± 100 

208 

3047 

2969 

^2960 ± 65 

3.6 GeV/ nucleon 

12 

836 

779 

**780 ± 30 

20 

1059 

902 

**1040 ± 60 

64 

1723 

1653 

**1700 i 90 


^Reference 17. 
^Reference 18. 
^Reference 19. 
^Reference 20. 


TABLE IV.- PAULI CORRELATION EFFECTS ON ABSORPTION CROSS 
SECTIONS FOR SELECTED COLLISION PAIRS 


Collision 

pair 

Incident 

kinetic 

energy^ 

GeV/nucleon 

^abs 

, mb 

Reduction 
in ^abs' 
percent 

Without 

correlation 

With 

correlation 

Ne + C 

3.6 

1144 

1059 

7.4 

0 + Cu 

.1 

2059 

1951 

5.2 


2.1 

2055 

1952 

5.0 


22.5 

2039 

1940 

4.8 

0 + Pb 

.1 

3416 

3341 

2.2 


2.1 

3418 

3346 

2.1 


22.5 

3403 

3333 

2.0 



TABLE V.- TOTAL CROSS SECTIONS FOR SCATTERING 



TABLE VI.- OPTICAL MODEL ABRASION CROSS SECTIONS 
FOR THE REACTION + 12q ^ n + X 

[Incident kinetic energy is 2.1 GeV/nucleon] 










Figure 5.*- Harmonic well charge and matter density 
distributions for 



r. fm 

Figure 6.- Harmonic well charge and matter density 
distributions for 



I 



r, fm 

• 12 

Figure 7*- C harmonic well charge distribution and the Woods-Saxon 
density obtained by artificially flattening it. 
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Figure 8.- Absorption cross sections for O projectiles. 
Incident kinetic energy is 2 .1 GeV/nucleon. 
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Figure 9.- Theoretical abrasion results (with and 

without the Pauli correlation correction) com- 

9 n 

pared with experiment for ^^Ne projectiles 
colliding with targets. Incident kinetic 
energy is 2 • 1 GeV/nucleon. 
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Figure 10.- Abrasion results for Ne projectiles 
colliding with targets predicted in this 

work compared with a previous theoretical model 
(ref. 10) and with experiment. Incident kinetic 
energy is 2 .1 GeV/nucleon. 
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